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ABSTRACT 

In an effort to understand the correlation between X-ray emission and present star 
formation rate (SFR), we obtained XMM- Newton data to estimate the X-ray luminosi- 
ties of a sample of actively starforming HII galaxies. The obtained X-ray luminosities 
are compared to other well known tracers of star formation activity such as the far 
infrared and the ultraviolet luminosities. We also compare the obtained results with 
empirical laws from the literature and with recently published analysis applying syn- 
thesis models. We use the time delay between the formation of the stellar cluster and 
that of the first X-ray binaries, in order to put limits on the age of a given stellar 
burst. We conclude that the generation of soft X-rays, as well as the Ha or infrared 
luminosities is instantaneous. The relation between the observed radio and hard X-ray 
luminosities, on the other hand, points to the existence of a time delay between the 
formation of the stellar cluster and the explosion of the first massive stars and the 
consequent formation of supernova remnants and high mass X-ray binaries (HMXB) 
which originate the radio and hard X-ray fluxes respectively. When comparing hard X- 
rays with a star formation indicator that traces the first million years of evolution (e.g. 
Ha luminosities) we found a deficit in the expected X-ray luminosity. This deficit is 
not found when the X-ray luminosities are compared with infrared luminosities, a star 
formation tracer that represents an average over the last 10^ years. The results support 
the hypothesis that hard X-rays are originated in X-ray binaries which, as supernova 
remnants, have a formation time delay of a few mega years after the starforming burst. 

Key words: Galaxies: evolution - Galaxies: starburst - X-rays: galaxies 



1 INTRODUCTION 

X-ray emission in star forming galaxies is dominated by a 
combination of high mass X-ray binaries (HMXBs) , hot O- 
stars, young supernova remnants and hot plasma, all sources 
that are closely related to the presence of massive short- 
lived stars which directly trace the current star formation 
activity (Persic & Rephaeli 2002) . Based on observations of 
local galaxies and the comparison with other tracers of re- 
cent star formation activity, several authors have proposed 
the use of the X-ray luminosity as a star formation tracer 
[e.g. David, Jones, & Forman (1992), Grimm, Gilfanov & 



* Partially based on observations obtained with XMM-Newton, 
an ESA science mission with instruments and contributions di- 
rectly funded by ESA Member States and NASA, 
t Visiting Fellow, loA, Cambridge, UK 



Sunyaev (2003), Ranalh, Comastri & Setti (2003)]. The lo- 
cally observed correlations appear to hold also at high-z and 
an empirical Lx — SFR relation based on observations of Ly- 
man Break Galaxies (LBGs) in the 1 Ms CDF-N has been 
derived (Nandra et al. 2002). However studies of LBGs suf- 
fer from the lack of other direct SFR tracers (e.g. Ha, FIR), 
and additionally, LBGs whose emission is dominated by star- 
burst events are too faint to be directly detected. Only the 
mean X-ray luminosities are attainable, making the possible 
contamination by obscured AGN difficult to quantify (Per- 
sic et al. 2004). Recent papers have confirmed the use of the 
X-ray emission as a direct tracer of star formation events 
[e.g. Laird et al. (2005), Rosa-Gonzalez et al. (2007a)]; in 
contrast, Barger, Cowie & Wang (2007) found a poor cor- 
relation between the X-ray and radio emissions of starform- 
ing galaxies that could question the use of X-rays as a re- 
liable tracer of star formation activity. In a recent paper 
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by Mas-Hesse, Oti-Floranes & Cervino (2008), the diflter- 
ent empirical relations between X-ray luminosities and the 
current SFRs were compared against evolutionary models. 
They found that under realistic assumptions (e.g. a young 
burst and an efficiency of a few percent in the rc-procossing 
of mechanical energy) the observed relations confirm the use 
of the soft X-ray luminosity as a reliable tracer of the star 
formation activity in young systems. 

This paper discusses the X-ray emission (as detected by 
XMM-Newton) of HII galaxies selected from the Tcrlcvicli 
et al. (1991) catalog as having intense Ha and If/3 emission 
lines in their optical spectra. 

HII galaxies are compact systems dominated by a strong 
and recent star formation burst event. The relative low mass 
- which implies a low contamination by low mass X-ray bina- 
ries - together with the absence of AGN activity make these 
objects the best laboratories to study the relation between 
the X-ray luminosity (Lx) and the current star formation 
rate (SFR). By using other tracers of star forming activity 
we will discuss the validity of the existing X-ray calibra- 
tions in very young systems. In fact, HII galaxies which are 
the youngest starbursts known in the local universe (Rosa- 
Gonzalez et al. (2007b) could have a deficit of X-ray emission 
due to a time lag between the formation of the massive star 
cluster and the formation of the first HMXBs. In the fol- 
lowing section we discuss the sample selection. In section 3, 
the data and its analysis. A comparison between the X-ray 
results and those obtained with other tracers of star forma- 
tion is presented in Section 4. Section 5 is the discussion and 
conclusions are given in section 6. 

Throughout this work a standard, flat ACDM cosmol- 
ogy with Q,A= 0-7 and Ho = 70 km s~^ Mpc~^ is assumed. 



2 SAMPLE SELECTION AND SFRS FROM 
OPTICAL EMISSION LINES 

HII galaxies harbour intense and young star formation 
events revealed by the strong emission lines observed in 
the optical spectra. Prom the Spectrophotometric Catalog 
of HII galaxies (Terlevich et al. (1991) we have selected 
those galaxies for which the corresponding SFR is greater 
than 4 M0yr~^. Based on the calibration given by Grimm 
et al. (2003) and the estimated SFR, we expect to find at 
least 12 HMXBs with luminosities greater than 10^* ergs~^ 
in each one of the selected galaxies. Due to the galaxies low 
mass, the expected number of LMXBs is a minimum. The 
selected sample is given in Table 1. 

The luminosities of the hydrogen recombination lines 
are proportional to the number of ionizing photons so they 
trace the presence of massive stars with lifetimes not larger 
than a few million years. We estimate the star formation rate 
from the reddening corrected Ha luminosities [SFR(Ha)] 
using the relation (Kennicutt 1998), 

SFR{Ha){MQyr-^) = 7.9 x W'^'^ L{Ha)ierg s'^) (1) 

Ha and H/3 fluxes are corrected for extinction using the 
Milky Way extinction curve (Seaton 1979) and assuming an 
intrinsic ratio Ha/H/3 of 2.87. The ionizing massive young 
clusters show in their spectra strong absorptions in the per- 
mitted lines (e.g. the HI Balmer series). The same species 
produces (at the same wavelengths) strong emission lines 



in the ionized gas and therefore the observed intensity of 
such lines is reduced by the underlying absorption, affecting 
the derived extinction. For those galaxies for which the H7 
line flux wEis available, we estimate the correction due to 
the presence of the low mass underlying population follow- 
ing Rosa-Gonzalez, Terlevich & Terlevich (2002). The final 
values of Av (the extinction in the V wavelength given in 
magnitudes) together with the Ha and H/3 fluxes and the 
derived SFRs are presented in Table 1. 

Notice that, as mentioned before, the SFR given by the 
Ha luminosity is tracing only the most massive stars (a few 
million years old), when a fraction of the expected HMXBs 
has yet to be formed. The time delay between the formation 
of the stellar burst and the formation of the first HMXBs, 
is one of the main topics of this paper and we will discuss it 
in detail in Section 4. 



3 XMM-Newton DATA ANALYSIS 

Fourteen objects have been observed with XMM-Newton 
(Jansen et al. 2001). The details of the observations are sum- 
marized in Table 2. The raw data have been processed using 
the standard Science Analysis System, SASv7.0.0 (Gabriel 
ct al. 2004). The most up-to-date calibration files available 
in July 2007 have been used for the data reduction. The 
raw data have been filtered from high background flaring 
using the method described in Piconcelli et al. (2004). De- 
tailed spatial and spectral analysis was possible only for four 
objects; net count rates were derived for all targets in the 
sample (Table 2). None of the objects was detected with the 
Reflection Grating Spectrometer. 

3.1 Imaging Analysis 

The majority of the sources are very weak in the X-ray 
band, allowing spatial analysis for only four of the galax- 
ies: Mrk 52, Tol 1457, Tol 1247, and Cam 0902. We have 
generated smoothed images for these galaxies by applying 
the SAS task asmooth^ to the pn 0.2-2 keV, 2-12 keV and 
0.2-12 keV band images (Figures 1, 2, 3, and 4). 

The images of Mrk 52 in different energy bands (Fig. 1) 
show a soft source, with no significant emission in the hard 
band. This is also the case of Tol 1247-232 (Fig. 4), although 
its emission in the hard band is stronger that in Mrk 52. The 
images of Cam 0902-1-1448 in the different energy bands are 
shown in Fig. 2. Tol 1457-262 is actually an HII galaxy pair. 
Both components are detected in the X-ray image but could 
not be resolved (Fig 3). The images show that the South- 
West member is the weakest one of the pair, with a stronger 
contribution in the hard band, i.e. 2-12 keV. 

3.2 Spectral analysis 

For the majority of the sources (ten out of fourteen) , the lim- 
ited signal-to-noise did not allow spectral analysis with the 
EPIC camera. However, enough counts were detected to pro- 
duce spectra for four sources: Mrk 52, Cam 0902, Tol 1457, 

■t The asmooth task was applied using the adaptive convolution 
technique with a S/N=8. 
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Table 1. Optical emission lines: fluxes, visual extinctions and SFHs. 



Name 


Distance 
(Mpc) 


F(Ha)(10-") 
(erg s~^cm~^) 


F(H/3) (10-") 
(erg s^^cm^ 


Av 


SFR(Ha) 
M0yr-i 


Mrk52 


31 


146. ± 22. 


17.0±2.5 


3.06±0.46 


14.42+4.84 


Cam0902+1448 


222 


28.1 ± 4.2 


3.40±0.51 


1.92±0.29 


50.48+10.92 


Toll457-262 


74 


85.8 ± 13. 


19.4±2.9 


1.50±0.22 


10.74+2.10 


Toll247-232 


213 


50.4 ± 7.6 


13.5±2.02 


0.74±0.11 


36.52+12.25 


Tol2259-398 


127 


19.9 ± 3.0 


3.09±0.46 


2.25±0.34 


14.95+5.02 


Cam08-82A 


222 


9.43 ± 1.4 


1 10+0 16 


3 05+0 461 


SS 45+1 2 90 


Mrk605 


131 


8.07 ± 1.2 


1.00±0.15 


2.87±0.43t 


10.15+3.40 


Tol2306-400 


292 


9.50 ± 1.4 


0.87±0.13 


3.71±0.56t 


106.99+35.89 


Mrk930 


78 


39.4 ± 5.9 


8.91±1.34 


1.20±0.18 


5.34+1.79 


UM530 


316 


26.9 ± 4.0 


4.37±0.65 


2.12+0.32 


114.40+38.37 


Tol0420-414 


87 


86.0 ± 13. 


28.2±4.23 


0.17±0.02t 


6.93+2.32 


T0IO6I9-392 


236 


9.53 ± 1.4 


1.23±0.18 


2.76+0.41 


35.72+11.98 


UM421 


177 


3.17 ± 0.48 


0.24±0.04 


4.24+0. 64t 


19.02+6.38 


UM444 


105 


17.1 ± 2.6 


3.24±0.49 


1.20+0.18 


4.15+1.39 



Av corrected by an underlying stellar population except when marked with a f 



Table 2. XMM-Newton observations ID, total and net exposure times after flaring removal, and background subtracted count rates in 
the diflferent X-ray bands. The data were extracted from the EPIC-pn camera. 



Target 


Obs. ID 


Exposure time 


Count Rate 




Count Rate 


Count Rate 










(0.15 - 10 keV) 


(0.5 - 2 keV) 


(2 - 10 keV) 






ks 




c/s 




c/s 


c/s 


Mrk 52 


0303560101 


6.84, 


4.3 


0.020+0.002 




0.0132+0.0018 


0.0008+0.0003 


Cam 0902+1448§ 


3035620101 


6.84, 


3.3 


0.110+0.011 




0.053+0.006 


0.030+0.008 


Tol 1457-262 


0303560601 


9.33, 


6.5 


0.059+0.004 




0.040+0.003 


0.013+0.003 


Tol 1247-232 


0303561001 


11.5, 


6.0 


0.024+0.003 




0.0150+0.0019 


0.0058+0.0014 


Tol 2259-398 


0303560701 


10.2, 


5.6 


0.010+0.003 




0.0057+0.0019 


0.004+0.003 


Cam 08-82A 


0303561101 


11.8, 


6.7 


<1.4xl0-3 




< 5 X 10-4 


< 6 X 10-* 


Mrk 605 


0303561701 


7.93, 


5.2 


(1.9+1.5)10- 


-3 


< 9 X 10-* 


< 9 X 10-5 


Tol 2306-400 


0303560401 


10.9, 


5.0 


(8+3)10-3 




(4.9+1.6)10-3 


0.003+0.002 


Mrk 930 


0303560901 


10.8, 


2.4 


(4.6+1.7)10- 


-3 


0.0017+0.0007 


0.0029+0.0009 


UM 530 


0303560501 


6.84, 


2.0 


(8+3)10-3 




(5+2)10-3 


< 4 X 10-4 


Tol 0420-414 


0303561901 


6.84, 


2.8 


(1.8+1.6)10- 


3 


(1.4+0.9)10-3 


< 3 X 10-4 


Tol 0619-392 


0303561401 


6.84, 


3.5 


(5.2+2.3)10- 


-3 


3.0+1.4)10-3 


0.0020+0.0015 


UM 421 


0303561601 


11.8, 


5.7 


(2.9+1.7)10- 


-3 


0.0017+0.0011 


0.0012+0.0009 


UM 444 


0303561801 


12.1, 


3.6 


(7+3)10-3 




(4.8+1.5)10-3 


< 8 X 10-4 



§ Another observation of Cam0902+1448 with ID 0303561201 is available. The results on spectral quantities derived from it, though 
less accurate, are compatible with those obtained from the 3035620101 observation so we used the latter for the analysis. 



and Tol 1247. Circulax extraction regions of 400, 700, 550 
and 650 pixels respectively were selected in order to maxi- 
mize signal-to-noiso in the 0.2-10 keV band (see details in 
PiconccUi ct al. (2004) . By using these regions we are includ- 
ing between 80% and 90% of the total energy as described in 
the XMM-Newton handbook. We also checked that we are 
not including X-ray emission from close sources and from 
pixels too close to the CCD edges. 

Background spectra have been extracted from a circu- 
lar region located in the frame of the galaxy and free of any 
visible contaminating source. The associated ancillary and 
response matrices were obtained using the standard SAS 
tasks. EPIC-pn and the combined MOS 1 and 2 spectra of 
the four brightest objects were extracted. The spectra of 
Cam 0902-1-1448, Tol 1457-262, and Tol 1247-232 have been 
binned such that each bin contains at least 20 counts in or- 
der to apply the minimization technique. As not enough 
counts were detected in Mrk 52, we used the unbinned spec- 
tra to perform the analysis and therefore C-statistics were 



applied. The XMM-Newton spatial resolution and the lim- 
ited signal-to-noise do not allow us to spectroscopically an- 
alyze each component of the pair Tol 1457 separately. Both 
members are included in the extraction region and therefore 
fluxes and luminosities refer to the pair's combined emis- 
sion. In all the galaxies, the spectral analysis has been per- 
formed using the vl2.0 of XSPEC. The errors quoted are 
referred to the 90% confidence level, i.e. Ax = 2.71 when 
statistics were applied. In all the cases, we fixed the 
Galactic column density to the values listed in Table 5 which 
were obtained through the Leiden/Argentine/Bonn Survey 
of Galactic HI (Kalberla et al. 2005). 

The high supernova rate occurring in star-forming re- 
gions produces an acceleration of the electrons to relativistic 
velocities. The scattering of these electrons with FIR pho- 
tons, enhanced by the star-forming process results in non- 
thermal X-ray emission. This emission can, at first order, be 
modelled as a power law with am index in the range 1.6-1.8. 
However, observational studies on star-forming galaxies find 
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Figure 1. Smoothed pn images in the 0.2-12 keV, 0.2-2 keV, and Figure 3. As in Figure 1, for Tol 1457-262. 

2-12 keV energy bands of Mrk 52. We also show the DSS image 
for comparison. Here and in figures 2, 3, and 4, North is up and 
East is left. 

EPIC-pn 0.2-12 keV EPIC-pn 0.2-2 keV 




Figure 4. As in Figure 1, for Tol 1247-232. 



Figure 2. As in Figure 1, for Cam 09024-1448. 

that a thermal emission model with temperatures of the or- 
der of 8 keV, provides also a satisfactory representation of 
the observed high energy spectra. The average X-ray spec- 
trum of a population of low- and high-mass X-ray binaries 
can be described by a power law with an index F ~ 1.2 and 
a cut-off around 7.5 keV (Persic & Rephaeli 2003). 

On the other hand, the diffuse emission mainly con- 
tributes to the soft X-ray band and can also be modelled 
as a thermal component but with a temperature one or- 
der of magnitude lower. Bearing this picture in mind, we 
have tested three different models to fit the source spec- 



tra: a single power law, a thermal emission model [specifi- 
cally the mekal model in XSPEC, Liedahl, Osterheld & Gold- 
stein (1995)], and a combination of the best of the two pre- 
vious ones with a second thermal emission model accounting 
for the soft emission. We also included a model for the photo 
electric absorption {zwabs in XSPEC) using the Wisconsin 
cross-sections (Morrison & McCammon 1983). 

The parameters fitted - by using a Levenberg-Marquart 
method - and the goodness of the fits are given in Table 3. 
The low number of detected counts does not allow us to fit 
more complex models. A single-component model provides 
acceptable fits to the data for all sources analyzed, except 
for Cam 0902, for which a two components fit is best. The 
derived fiuxes and unabsorbed luminosities for each object 
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are given in Table 4. In Figs 5 and 6, we illustrate the ob- 
served spectra and the best fit models for each gala^cy. To 
illustrate the variation of the errors in the T vs. Nh and 
the temperature vs. Nh planes we plotted in Figure 7 the 
corresponding error ellipses for Cam 0902. 

We compare the obtained X-ray fluxes (Table 4) with 
those given by a starburst model described below (Table 5). 
The agreement is good but some discrepancies were found 
in the soft X-ray fluxes mainly due to the uncertainties in 
the Nh value. In fact, differences in the goodness of the flt 
were found when the fit was made by fixing the value of the 
Nif derived from the Hydrogen emission lines {Nh = 1-79 x 
10'^"'^ X Av, cm"'^) as compared with the results obtained 
when the Nh is a free paranietcr of the fitting procedure. 
The derived values of F are somewhat higher than those 
found in the literature for similax type of objects and fix 
F = 2 created fits with large values of x^- However the 
data quality does not allow the discussion of the origin of 
this apparent systematic effect. In one case (Cam 0902), the 
best fit was obtained by adding to the power law component 
a thermal model with a temperature of about 0.6 keV. That 
is consistent with temperatures of the hot plasma observed 
in other star forming galaxies (Grimes et al. 2005). 

In order to have a homogeneous data set, and due to 
the low rmmber of photons that could be affecting the con- 
clusions obtained from the spectra, the analysis presented 
in the next sections is based only on the count rates from 
Table 2. The count rates were calculated for the pn detector, 
extracted from a region of 15" free of contaminating sources 
and corrected for background. We convert the count rates 
by using a two component model; a power law with spectral 
index of 1.2 and a cut-off at 7.5 keV (associated to the com- 
bined high and low-mass X-ray binaries contribution) plus a 
thermal emission model with kT=0.8 keV (associated to the 
diffuse emission from galactic winds). The two components 
(power law and thermal emission) arc normalized such that 
the fiux ratio between the thermal and the power law com- 
ponents in the 2 — 10 keV band is 0.03 [(Persic & Rephaeli 
2003), and references therein]. The estimated fiuxes in the 
soft (0.5-2 keV) and hard (2-10 keV) bands are shown in 
Table 5. Calculations were performed using P/MA/5'v3.9^. 

The fluxes and luminosities, also given in Table 5, were 
corrected from absorption using the extinction derived from 
published ratios between Ha and H/? which include the ab- 
sorption due to the Galaxy. Notice that the values of the 
cohnnu density given in Table 5 were used only in the spec- 
tral analysis of Mrk52, Cam0902, Toll457 and Toll247. The 
errors were derived from the count rate uncertainties. How- 
ever, we notice that they should be considered lower limits 
of the real uncertainties, as they only include the poisson 
error on the count rate. The uncertainty generated by con- 
sidering different models is not included in the error budget. 
In this sense, when a single power-law model with F = 2 
is applied instead of the cut-off power-law plus thermal 
emission, as suggested by other studies on single sources 
[e.g. NGC 3310, NGC 3690 by Zezas, Georgantopoulos, & 

§ PIMMS docs not consider the cut-off power-law for 
fluxes/count rates estimates. Therefore, we have also applied a 
correction due to differences on flux estimation when a simple 
power-law model is used instead of a cut-off power law. 



Ward (1998); NGC 253, M82 by Cappi et al. (1999)], differ- 
ences of the order of 20 percent in the estimation of fluxes 
are obtained. The metallicity of the gas is another source of 
uncertainty particularly in the soft band. Calculations real- 
ized with PIMMS show that for a plasma with a temperature 
of 0.8 keV a change in metallicity from solar to 20% solar, 
typical of HII galaxies [e.g. Hoyos & Diaz (2006)], produces 
about 15% variation in the fluxes. Finally, we estimate the 
luminosity due to low mass X-ray binaries (Llmxb) that 
could be present in the selected galaxies. Based on 2MASS 
K-band magnitudes we confirm that the observed galaxies 
have masses which are in general lower than that of an L* 
galaxy (assuming a mass to light ratio in the K-band of 0.5, 
Bell & dc Jong (2001). The contribution of these binaries 
to the observed X-ray luminosities is estimated using the 
relation between the stellar mass and the luminosity due to 
LMXB (Grimm, Gilfanov & Sunyaev (2002), 

Llmxb {erg s'^) = 5 x 10^* (2) 

where M« is the stellar mass in solar units. Table 5 shows 
that the calculated Llmxb are much smaller, on average less 
than 2 percent and at most 20 percent, than the observed 

X-ray luminosities. The general result is that the total X- 
ray luminosities of the low mass galEixies selected here is 
probably not greatly affected by the LMXBs contribution. 

4 EARLY EVOLUTIONARY PHASES: 
COMPARISON WITH SFR TRACERS 

4.1 Hq luminosities 

The recombination emission lines only appear during the 
early evolution of a starburst when the massive stars are hot 
enough to produce ionizing photons. Therefore, the deduced 
SFR(Ha) is tracing the first million years of the present star- 
burst evolution. Figure 8 shows the relation between the Ha 
luminosity and the X-ray luminosities both in the soft (top 
panel) and in the hard (bottom panel) bands. The solid lines 
show the empirical relation found by Ranalli et al. (2003) 
based on the study of nearby star forming galaxies, 

SFRsoft{MQyr-^) = 2.2 x 10"''°L(0.5-2fecy) (erps"^)(3) 

SFRhard{MQyr-^) = 2. x 10"''"L(2 - lOfceV) (erps"')(4) 

There is a good correlation between the L(Hq) and the 
soft X-ray luminosity (top panel), but there is one galaxy 
(Cam 0902-1-1448) for which the luminosity in X-rays is well 
above the one predicted by Ranalli's law. In the bottom plot 
of this Figure [hard X-rays vs. L(Hq)] most of the galax- 
ies are below the empirical correlation. Due to the pres- 
ence of upper limits we fit the luminosities by using survival 
analysis^. The Buckley- James algorithm [BJ, e.g. Isobe, 
Feigelson, & Nelson (1986)] a non-parametric method which 
treats the residuals using the Kaplan-Meier description, was 
applied to the data in Figure 8. In the case of the bot- 
tom panel, the non-parametric method was applied both to 
all the data (dashed — line) and excluding Mrk52 (dotted — 
line) . 

^ See IRAF task stsdas.analysis.statistics. survival for an intro- 
duction to survival analysis. 
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Figure 5. Observed EPIC spectra, best fit model and residuals of left panel: Mrk 52; right panel: Cam 0902+1448. 
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Figure 6. As in Figure 5 for left panel: Toll247-232; right panel: Toll457-262. 



4.2 Radio luminosity 

Table 6 shows the radio fluxes at 1.4 GHz for the selected 
galaxies extracted from NED. The strong correlation known 
to exist between the radio and the far infrared fluxes is taken 
as an indication that radio emission is a reliable tracer of 
the SFR. Yun, Reddy, & Condon (2001) found, from a FIR 
selected sample of about 1800 star-forming galaxies, that 
their radio continuum may be used to infer the extinction- 
free SFR using the relation 

SFR{lAGHz){MQyr'^) = 5.9x10'^^ L{lAGHz){erg s'^ Hz' 

The values of SFR(1.4GHz) for our XMM-Newton observed 
galaxies are also presented in Table 6. The empirical rela- 



tion found by Yun et al. (2001) is well understood and it 
is based on the existence of a star forming event in which 
the observed radio emission arises from accelerated electrons 
produced in the supernova remnants and the IR radiation 
comes from dust that is heated mainly by massive stars. The 
estimated SFR corresponds to continuous star formation oc- 
curring for the last 10* years. 

In Figure 9 we compare the radio with the soft X-ray 
luminosities together with the empirical relation given by 
Ranalli et al. (2003). In the left panel the X-ray luminosities 
)|^e calculated by using a two component model (thermal -|- 
power law as described in section 3.2). Although the original 
calibration by Ranalli et al. (2003) is based on the relation 
between the radio and X-ray luminosities (solid line) , a clear 
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Table 3. Spectral analysis. 



Target 


Model 




r 




kT 


Goodness 




(1) 
















lO^^^cm^ 




10^1 cm^ 


keV 




Mrk 52 


A 


6 ± 3 


6.5 it 0.5 






489 tor 3222 dot+ 




At 


5.5f 


6.5 ± 0.5 


; 


; 


490 for 3223 dof* 




B 


- 


- 


7 7+1.5 
'•'-1.3 


0.14 ±0.02 


497 for 3222 dof* 




Bf 






5.5f 


1O+0.02 
U-lS-0.03 


505 for 3223 doft 


Cam 0902+1448 


A 


^■°-1.2 


c 9+0.8 


- 


- 


31 for 25 dof 




At 


3.4f 


r, O+0.3 


- 


- 


35 for 26 dof 




B 






0.8 


0.2 


DO tor 25 dot 




Bf 






3.4f 


0.25 


73 for 26 dof 




C 


2 0+^-^ 
^•'^-1.6 






63+0 i« 


18 for 23 dof 


Toll457-262 


A 


1 4+1-0 


r, 9 + 0.5 




4 6+3 


16.2 for 12 dof 




B 






< 0.6 


18.7 for 12 dof 


Toll247-232 


A 


< 1.4 


1 S+0'4 
^■^-0.3 






5.1 for 6 dof 




B 




< 0.5 


7+16 
'-4 


7.6 for 6 dof 



Notes: (1) Model A: zwabs*zpowerlaw;Modcl B; zwabs*mckal; Model C: zwabs(zpowerlaw+mckal) ; in boldface, the best fit model. | Nh fixed to the value 

derived from Ay. ^ C-statistics. 



Table 4. Observed fluxes and unabsorbed luminosities in the soft and hard bands. 



Target Flux (0.5-2 keV) Flux (2-10 keV) Luminosity (0.5-2 keV) Luminosity (2-10 keV) 

10-11 erg s-i lO^'i erg s^i 



Mrk 52 2.5 ±0.5 

Cam0902+1448 8.5^4^ 



Toll457-262 
Toll247-232 



5.5 
2.2 



+ 1.3 

1.9 
+0.8 
1.3 



0.08±0.03 

1-0-1.3 



5.8 ±1.2 

1+2 
0.7 ±0.2 

K+0.6 

-0.9 



14"' 



1.6"' 



1.0 ±0.4 X 10^ 




0.2 0.4 
NH(xlO==cm-=) 




0.2 0.4 

NH(xlO==cm-=) 

Figure 7. Confidence contours for the determination of F and 
A'^j^ (top panel) and temperature and Nh (bottom panel) for 
Cam0902. In both panels the best fit is marked with a cross. 



10" 
10" 
10" 
10^ 

10" 
10" 
10" 
10^ 




Cam0902 



'MrkoZ 



Cam0902 




10" 



10" 



L„„ (org s"') 



Figure 8. X-ray (soft band in the top and hard band in the bot- 
tom panel) vs. Ha luminosities. Here and in the following plots, 
solid circles mark those galaxies for which we extracted the X-ray 
spectrum. In both cases, the solid line is the result of converting 
the Ranalli law to Ha luminosities and the dashed one is a linear 
fitting that takes into account the presence of upper limits in the 
derived X-ray luminosities. In the bottom panel, the dotted line is 
the result of the fit obtained when leaving out the outlier Mrk52. 
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Table 5. Galactic column density and derived fluxes and luminosities. 



Target Nj/ Fo,5_2fcel' F2-10fceV' L().5_2fceV L2_10feeV ^l)^^™keV ^LMXB 

(lO^O) (10-14) (10-14) (1041) (1041) (1041) (1039) 

cm"^ erg/s/cm erg/s/cm erg/s erg/s erg/s erg/s 



ivirK oz 


1 

i, 


o 
.y 




y.yUiti.U4 


U.lUit U.Ui 


n 1 1 -1- n ni 
U.liit U.Ui 


U.U4oitU.UU0 


1 

1, 


. / D±U. 


OO 


Cam 0902+1448 


,§ 3, 


.7 


41.00±4.13 


48.00±4.84 


24.24± 2.44 


28.38± 2.86 


12.1 ±1.22 


2, 


.34±0, 


,70 


Tol 1457-262 


9, 


.4 


20.00±1.36 


24.00±1.63 


1.30± 0.09 


1.56± 0.11 


0.651±0.044 


0, 


.47±0, 


14 


Tol 1247-232 


6, 


.6 


6.40±0.77 


7.50it0.90 


3.48± 0.42 


4.07± 0.49 


1.73±0.209 




< 0. 


,73 


Tol 2259-398 


1, 


.1 


3.70±1.11 


4.30±1.29 


0.71± 0.21 


0.83± 0.25 


0.357±0.107 


1, 


.04±0. 


,31 


Cam 08-82A 


3, 


.3 


<0.43 


<0.50 


<0.25 


< 0.30 


< 0.127 




< 0. 


,63 


Mrk 605 


1, 


.9 


<0.82 


<0.97 


<0.17 


< 0.20 


< 0.085 


1, 


.36±0. 


,41 


Tol 2306-400 


1, 


.3 


3.73±1.40 


4.30±1.61 


3.81± 1.43 


4.39± 1.65 


1.90 ±0.714 


3, 


.96±1. 


,19 


Mrk 930 


5, 


.5 


1.40±0.52 


1.61±0.60 


O.lOit 0.04 


0.12± 0.04 


0.051±0.019 




< 0. 


,08 


UM 530 


1, 


.5 


2.74±1.03 


3.18±1.19 


3.27±1.22 


< 3.79 


1.63 ±0.612 


6, 


.43±1, 


,93 


Tol 0420-414 


2, 


.5 


0.58±0.52 


0.69±0.61 


0.05±0.05 


< 0.06 


0.026±0.023 




< 0, 


13 


Tol 0619-392 


7, 


.6 


2.20±0.97 


2.60±1.15 


1.47± 0.65 


1.73± 0.77 


0.734±0.325 


1, 


.78±0, 


,53 


UM 421 


4, 


.1 


1.50±0.88 


1.80±1.06 


0.56± 0.33 


0.67± 0.39 


0.280±0.164 


2, 


.28±0. 


,69 


UM 444 


2, 


.1 


2.20±0.94 


2.60±1.11 


0.29±0.12 


< 0.34 


0.144±0.062 


0, 


.54±0. 


,16 



deviation from the empirical relation is obtained for our ob- 
jects. Most of the gala^xies are below the line defined by 
Ranalli and collaborators and seem to have low radio emis- 
sion. These galaxies could be in an early stage previous to 
the synchrotron phase that appear with the first supernova 
explosions. At this phase the radio emission is dominated by 
free-free electrons produced in HII regions close to massive 
stars. 

To check if a better correlation is found when only the 
X-ray thermal emission is considered we use a thermal model 
to convert the counts to luminosities. Wo used PIMMS and 
a thermal mekal kT=0.8 keV model with solax metallicity 
to derive the luminosities for all the galaxies in the sample. 
In this model we assume that wo are counting photons pro- 
duced by the interaction between the hot out-flowing wind 
and the ambient gas in the host galaxy (Grimes et al. 2005) . 
In the right panel the X-ray luminosities correspond to the 
so defined thermal component only. In both panels the solid 
line represents the empirical relation of Ranalli and collabo- 
rators and the dashed hne, the linear regression to our data. 

The radio fluxes of flve of the galaxies axe just upper 
limits, so we implement the BJ algorithm to confirm the 
observed trends. The results axe shown by the dot-dashed 
line in both panels of figure 9 and in figure 10. 



4.3 IR and UV luminosities 

As wo mentioned before, the IR radiation is a good tracer 
of the SFR averaged over the last 10* years. However an 
unknown fraction of UV photons could have escaped the 
galaxy and consequently, the derived SFR would be under- 
estimated. To solve this problem Heckman et al. (1998) and 
Buat et al. (1999) proposed to add the SFR calculated from 
the ultraviolet to that based on the IR. By adding the two 
contributions, on one hand one avoids the complex correc- 
tions that must be applied to the UV data and on the other 
hand, there is no need to correct for the uncertain UV es- 
caping fraction that affects the IR measurements. We use 
the relations (Kennicutt 1998), 

SFR{IR){MQyr~'^) = 4.5 x 10"^'^ L{IR){ergs~^) (6) 



SFR{UV){MQyr~'^) = 1.4 x W-'^^L^{ergs~^Hz~^) (7) 

which are valid for the case of a universal Salpeter initial 
mciss function with masses between 0.1 and 100 M©. The 
UV data comes from GALEX (MAST Archive) and our UV 
observations made with OM on board XMM-Newton (Ta- 
ble 6). For those galaxies for which IRAS data exist, we 
estimate both the SFR(IR) and SFR(UV). The total SFR 
for them is calculated directly by adding the SFR(IR) and 
the SFR(UV) without applying any extinction correction. 
The results arc given in Table 6. For those galaxies lacking 
IRAS data, the total SFR is calculated from the UV fluxes 
corrected by extinction, using Equation 7. To correct for ex- 
tinction we use Calzetti's law (Calzetti 2001) and the visual 
extinctions calculated from the Balmer decrement (Table 1) 
but taking into account that due to gas surrounding the 
newly formed stars, they are more extinguished than the 
stars from which most of the UV radiation is coming [age 
selective extinction, e.g. Mayya et al. (2004)]. Cid Fernan- 
des et al. (2005) based on the study of 50,362 galaxies from 
the Sloan Digital Sky Survey (SDSS) found that the relation 
between the visual extinction obtained from the Hq to H/3 
ratio is Av = 0.24 -|- l.SlAv* , where Av* is the extinction 
that is affecting the continuum light. The Av* values axe the 
ones that we apply to correct the observed UV fluxes. The 
comparison of the total SFR with the X-ray luminosities is 
presented in Figure 11. The ratio between the soft X-ray 
and IR luminosities can be used to study the evolutionary 
stage of a given star forming burst. The values are given 
in Table 6 and a comparison with the results from recent 
theoretical models is discussed in the next section. 



5 DISCUSSION 

The strong hydrogen recombination lines observed in the 
selected galaxies suggest the presence of an intense burst 
younger than a few million years. In most of the cases, the 
derived SFR is higher than that of the archetypal starburst 
galaxy M82 (Mayya et al. (2004). In the early phases of 
the evolution of a star cluster (just before the flrst SN ex- 
plosions) the soft X-rays are mainly produced by the me- 
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Figure 9. Radio versus soft X-ray luminosity (Lg.s— 2fceu) for the selected galaxies. The X-ray luminosities have been calculated by 

using a two model component (thermal + power law) in the left panel and a thermal component in the right panel. Galaxies for which 
it was possible to extract the spectra arc marked with solid circles, the other sources in our sample are indicated as crosses. In both 
panels, the galaxies studied by Ranalli et al. are represented by asterisks. The solid line shows the linear correlation found by Ranalli and 
collaborators. The result of a standard linear fit (dashed line), and the fit from the survival analysis (dot — dashed line) are also plotted. 



Table 6. Fluxes and SFRs from radio, infrared, ultraviolet and total. The ratio between the soft X-ray and IR luminosities is shown in 
the last column. 



Name 


F(1.4GHz) 


SFR(1.4GHz) 


F(60 /im) 


F(100 i^m) 


SFR(IR) 


F(UV) 


SFR(UV) 


Tot. SFR 


1°S L(IR) 




(mjy) 


Moyr-i 


(Jy) 




(Mgyr-i) 




(Moyr-i) 


(Mgyr-i) 


Mrk52 


13.10±0.60 


1.10±0.33 


4.43±0.03 


6.65±0.10 


1.47± 0.44 


3086.38±36 


0.61± 0.31 


2.09± 1.2 


-3.51 


Cam0902-|-1448 


86.84±0.82 


309±93 


4.12±0.21 


6.98±0.35 


62.1± 19. 


1207.58±24 


10.41± 5.2 


72.6± 42 


-2.76 


Toll457-262 


37.80±1.80 


14.5±4.4 


3.09±0.18 


3.68±0.41 


4.34± 1.3 


1670.91±21 


1.52± 0.76 


5.87± 3.4 


-2.87 


Toll247-232 


3.40±0.50 


10.9±3.6 


0.51±0.05 


< 0.97 


10.6± 3.2 


943.07±12 


7.17± 3.6 


17.8± 10. 


-2.83 


Tol2259-398 


< 1.77 


2.01 








477.96±3.4 


1.29± 0.64 


lO.lit 5.0 




Cam08-82A 


3.06±0.15 


10.7±3.2 


0.58±0.01 


0.95±0.16 


8.32± 2.5 


203.48±3.9 


1.68± 0.84 


10. 0± 5.8 




Mrk605 


< 1.56 


1.90±4.3 


0.27±0.04 


< 1.36 


2.17± 0.65 


172.18±4.3 


0.50it 0.25 


2.67± 1.6 




Tol2306-400 


< 1.38 


8.31 


<0.25 


< 0.90 


<9.86 


325.57±11 


4.65± 2.3 


162± 81 




Mrk930 


12.20±0.90 


5.26±1.6 


1.25± 0.09 


< 2.15 


3.52± 1.1 


864.69±18 


0.88± 0.44 


4.40± 2.6 


-3.89 


UM530 


6.40±0.60 


45.0±14 


0.58±0.07 


0.63±0.13 


14.5± 4.3 


340.44itl7 


5.68it 2.8 


20.1± 12. 


-2.99 


Tol0420-414 


< 1.50 


0.80 


<0.25 


< 0.90 


<0.87 


160.44±5.6 


0.20± 0.10 


0.20± 0.1 




T0IO6I9-392 


< 1.56 


6.14 


0.42±0.05 


< 1.59 


10.8± 3.2 


238.70±7.8 


2.23± 1.1 


13.0± 7.6 


-3.21 


UM421 


4.30±0.50 


9.46±3.0 


0.56±0.07 


< 1.30 


8.07± 2.4 


101.48±7.7 


0.53± 0.26 


8.60± 5.0 


-3.51 


UM444 


< 1.44 


1.11 


<0.25 


< 0.90 


<1.27 


465.66±7.8 


0.85± 0.42 


2.28± 1.1 





chanical energy released by stellar winds [e.g. Strickland 
et al. (2002), Corvifio, Mas-Hesse, & Kunth (2002), Silich, 
Tenorio-Taglc, & Afiorvc-Zofcrino (2005)], when the Ha lu- 
minosity is also stronger (Figure 8 top panel). The hard 
X-ray luminosity, on the other hand, is generally attributed 
to HMXBs (e.g. Grimm et al. (2003) which, at the early 
stages of evolution of these clusters are still in the forma- 
tion process. In fact the maximum number of HMXBs occurs 
20-50 Myr after the star forming event (Shtykovskiy & Gil- 
fanov (2007). Therefore, a lack of hard X-rays is expected 
with respect to the instantaneous star formation rate traced 
by the Ha luminosity (Figure 8 bottom panel). 

A dearth of radio emission in the early phases of clus- 
ter formation has been attributed either to a strong free-free 
absorption of the emitted radio waves (more severe at longer 
wavelengths), or to an epoch when most of the massive 
stars have not yet exploded as supernovae [Rosa-Gonzalez at 
al. (2007b)]. Ha and hard X-ray luminosities probably dif- 
fer because of the time delay between the formation of the 



first SN remnants and the emergence of stellar winds in the 
cluster massive stars. Therefore it is hard to interpret SFRs 

based on different tracers which in general arc related to dif- 
ferent time scale phenomena. Figure 9 (radio vs. soft X-ray 
luminosities) shows the clear signature of a young starburst, 
where the radio luminosity is well below that expected from 
the Ranalli ct al. (2003) relation. We have found that the 
SFR(1.4GHz) is in general lower than the SFR(Hq) with 
the remarkable exception of Cam 0902-1-1448 for which the 
SFR(1.4GHz) is about 6 times larger than the SFR(Ha). 
When wo compare the radio luminosities with the soft X-ray 
luminosities and the Ranalli empirical relation the deficit in 
radio is noticeable (Figure 9). However, when we compare 
the radio luminosities with the fraction of the soft X-ray 
luminosities due to thermal emission alone, a better agree- 
ment is found (sec right panel of Figure 9). Notice that when 
the BJ method is used to estimate the linear regression, a 
higher discrepancy is observed. 

Radio versus hard X-ray luminosities are shown in Fig- 
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urc 10. We can see that most of the galaxies are below the 
relation proposed by Ranalli et al. (2003), consistent with 
the lack of both SN remnants and HMXBs. 

For galaxies with available IRAS data, the total SFR 
was calculated by adding the SFR(IR) and the SFR(UV) 
without correcting the UV fluxes for extinction. For those 
galaxies without IRAS fluxes, the total SFR was calculated 
from the extinction corrected UV luminosities. The total 
SFR is an average over the last 10* years, providing enough 
time for a population of HMXB to be in place. In fact there is 
a good correlation between the total SFR (averaged over the 
last 10* years) and the X-ray luminosities both in the soft 
and in the hard band (Figure 11), in fairly good agreement 
with the empirical Ranalli relation. 

Notice that all selected galaxies are classified as HII 
galaxies based on BPT diagrams of optical line ratios. How- 
ever an unknown contribution from a heavily obscured AGN 
could be present. This contribution is more important in the 
hard band, and the deviation from the Ranalli relation ob- 
served in Figure 10 could be due to an excess of X-ray emis- 
sion from an AGN instead of representing a deficit in radio. 
However, such deficit in radio is clear when the SFR(Ha) 
is compared with the SFR(1.4 GHz). Moreover, the excess 
of hard X-ray emission is not present when comparing the 
hard band X-ray luminosity with the total SFR (Figure 11). 
We can safely conclude that the X-ray emission that we see 
is due to star formation processes. 

Mas-Hesse et al. (2008) have concluded that assuming 
a reprocessing efficiency of the mechanical energy of a few 
percent, the observed values of the ratio between the soft X- 
ray and infrared luminosities can be reproduced with syn- 
thetic models of young starbursts. In fact they obtain an 
average value of log Loj^-'zkeV / Lm — —3.5 for efficiencies 
between 1% and 10% for bursts with ages around 5 Myr. 
In our sample, the ratio log Lo.5-2kev /Lm (Table 6) has an 
average value of —3.20 ± 0.41 consistent with the theoret- 
ical estimations. For continuous star formation, a value of 
log I/o.5-2fccV' /LiR = —3.5 is consistent with ages as large 
as 25 Myr. At this age, a significant number of stars have 
exploded as SN. The strong emission lines and the high SFR 
derived together with the lack of synchrotron emission ob- 
served in most of the sources suggest a star formation history 
characterized by a strong short burst. 

Mrk 52 and Cam0902 are outliers in some of the pre- 
sented plots. Particularly, Mrk 52 is a bright blue galaxy 
with a derived SFR(Hq) = 14 Moyr"^. However, the ex- 
tremely weak emission in the hard X-ray band (see Figure 1) 
together with the deficit in radio makes this galaxy the most 
clear example of an extreme young burst with an age of 
about 3 Myrs when the majority of supernovae and HMXBs 
have not yet been formed. This age is consistent with Mrk 52 
being a Wolf-Rayet galaxy as the dominant burst is going 
through an intense Wolf-Rayet phase [Schaerer, Contini, & 
Pindao (1999), Fernandes et al. (2004)]. 



6 CONCLUSIONS 

We have analyzed XMM-Newton observations of a sample of 
14 star forming galaxies with high SFR judging from their 
strong hydrogen recombination lines. The main conclusions 
of our analysis are: 




Figure 10. Radio vs. hard X-ray luminosities. Symbols as in 
Figure 9. 




0.1 1.0 10.0 100.0 

Total SFR 



Figure 11. X-ray luminosity (top panel for the soft and bottom 
panel for the hard band) against the total SFR. The line shows 
the Ranalli empirical law for soft and hard X-ray luminosities. 
In both panels, solid circles are the bright X-ray objects in the 
sample with an X-ray spectrum. 



(i) The SFR obtained from the soft X-ray luminosity is 
comparable to that determined from the Ha luminosities. 
This is consistent with the X-ray luminosity as a good tracer 
of SFR even at early times when the massive stars dominate 
the energy generation. Due to the time delay between the 
formation of the burst and the formation of the first massive 
binary systems, a lack of hard X-ray luminosity is observed 
when compared with tracers of recent star formation activity 
(e.g. Ha, Figure 8). This effect has already been observed in 
the Small Magellanic Cloud; we have extended the study to 
a larger sample of starforming galaxies. 

(ii) The weak radio emission observed in some of the ob- 
jects (related to early phases in the starburst evolution) to- 
gether with the values found for log Lo.5-2keV /Lir ~ -3.2 
suggests that the sample of galaxies is biased to galaxies 
dominated by young bursts with ages lower than 5 Myrs 
and with star forming histories characterized by long and 
low activity periods, followed by strong short episodes. 

(iii) The relation of both soft and hard X-ray luminosity 



© 2009 RAS, MNRAS 000, 1-11 



Evolution of the X-ray Luminosity in Young HII Galaxies 11 



with the SFR traced by the IR and UV, shows that the 
X-ray production is maintained during at least 10* years. 
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